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ABSTRACT 
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Q>^ , MHD turbulence is generally believed to have two important functions in accretion 

disks: it transports angular momentum outward, and the energy in its shortest 
wavelength modes is dissipated into the heat that the disks radiate. In this paper 
we examine a pair of mechanisms which may play an important role in regulating 
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the amplitude and spectrum of this turbulence: photon diffusion and viscosity. 
We demonstrate that in radiation pressure-dominated disks, photon damping of 
compressive MHD waves is so rapid that it likely dominates all other dissipation 
mechanisms. 
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1. Introduction 



Turbulence is widely thought to be central to the dynamics of accretion disks. A combination 
^ I of magnetic and Reynolds turbulent stresses may be responsible for the outward transport of 

^ ' angular momentum without which no accretion could occur ( ^hakura and Sunyaev 1973 , Balbus 



et al. 1994). The energy put into this turbulence is ultimately deposited as heat, and is therefore 
the energy source for the radiation by which we observe accretion disks. Although much effort 
has gone into identifying mechanisms which excite turbulence ( [Balbus and Hawley 199l| ), far less 
attention in the literature has been given to how the turbulence dissipates. In most instances, it is 
simply assumed that nonlinear couplings transfer energy from long wavelengths to short, and that 
some dissipative mechanism eventually damps very short wavelength motions. 

One reason why little thought has been given to the specifics of dissipation is that, as matter 
drifts inward through an accretion disk, if the disk is in a time-steady state its lost potential 
energy is transformed into heat and kinetic energy at a rate which is entirely fixed by global 
properties. If the gravitational potential is dominated by the mass M of the central object, the 
heating rate per unit area is 

3 GMM ^ , ^ 

Q = ^^:^RRir). (1) 

Rr (~ 1 at large radii) describes the reduction of the local heating due both to the kinetic energy 
carried outward with the angular momentum flux, and relativistic effects should the central object 
be a neutron star or black hole ( |Novikov and Thorne 1973 ). 
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It is a great simplification to calculations of disk equilibria that the heating rate should depend 
only on global quantities. However, this fact leaves open the question of how exactly the energy 
lost by the accretion flow is transformed into heat, and there are strong observational consequences 
that depend on just how this happens. For example, the existence of weakly-radiative disks 
( Ichimaru 1977] ; Rees et al. 198^ ; Narayan and Yi 1995) depends critically on the assumption 



that most of the heat goes to the ions, not the electrons. There have been other suggestions that 



a significant part of the heat goes into non-thermal particle distributions (e.g. [Ferrari 1984| or 



Stecker et al. 1991 ). Alternatively, the energy can be lost in magnetic fields which escape the disk. 



forming a corona or outflow (Galeev et al. 1979). 

Balbus & Hawley (1991) pointed out that MHD fluctuations should be linearly unstable 
in weakly-magnetized accretion disks. Fully nonlinear simulations (Hawley, Gammie & Balbus 
1995; Brandenburg et al. 1995; Stone et al. 1996) have shown that these fluctuations grow 
until the field energy density approaches the pressure in the disk, and that nonlinear couplings 
create fluctuations on shorter and shorter wavelengths. Most recent work on how the energy in 
these fluctuations is dissipated has concentrated on plasma physics effects that work on modes of 
very short wavelength (e.g. Bisnovatyi-Kogan & Lovelace 1997; Quataert 1997; Blackman 1997; 
Gruzinov 1997), especially in low-density, high temperature disks. 

Although this focus is well-grounded in reality in the context of MHD turbulence in laboratory 
plasmas, it ignores the fact that accretion disks are often extremely bright, and can contain such 
high photon densities that radiation dominates the total pressure. In this paper we point out 
that photon diffusion and viscosity can, in radiation- dominated accretion disks dominate all other 
mechanisms of dissipation. When that is so, compressive modes whose wavelengths are almost as 
great as a disk thickness can be rapidly damped. Significant consequences follow for the amplitude 
of MHD turbulence, the rate at which angular momentum may be transported, and the way in 
which the energy associated with the turbulence is dissipated into heat. 

The structure of this paper is as follows: we first extend (§2) the theory of MHD modes 
interacting with a background photon gas by substituting a time-dependent radiation transfer 
solution for the conventional description in terms of a photon viscosity. Our procedure is similar 
in character to the one adopted to treat photon diffusion damping of perturbations in the early 
Universe ("Silk damping": Silk 1968, Hu & Sugiyama 1996). We then apply this improved theory 
to conventional accretion disk models (§3). In §4 we discuss the impact of photon damping on 
both advection-dominated accretion disks and disks in which the dissipation is segregated into a 
corona. Finally, in §5 we summarize our results and discuss their significance. 

We close this introduction with some notes of distinction. There were earlier suggestions by 



Loeb and Laor (1992) and Tsuribe and Umemura (1997) that photon viscosity due to an external 



radiation field might explain the radial angular momentum transport in some accretion disks. We 
do not make that claim; in this paper we consider only how photon kinetic effects help regulate the 
amplitude of the MHD turbulence that is responsible for angular momentum transport. The effects 
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of photon damping we consider are for a scattering-dominated plasma, and thus the relations we 
derive are different from those found by Bogdan and Knolker (1989) and Mihalas and Mihalas 
(1983), who derived the dispersion relation for a radiation field in LTE, ignoring both scattering 



and radiation viscosity, which we include. Our problem also differs from that treated by [Cassen 
and Woolum (1996), who considered only optically thick spiral density waves that lose angular 



momentum through radiation. Our equations are very similar to those of Jedamzik et al. (1998) 



and Subramanian and Barrow (1997) in the diffusion and free-streaming limits; however, we have 
bridged the two regimes by truncating the radiation field moment expansion above quadrupole 
moment. We also note that Thompson and Blaes (1998) have considered radiation damping for 
waves in the context of gamma ray bursts. 



2. Equations 

Our aim in this section is to derive a dispersion relation for MHD waves in the presence 
of a background radiation field. In a sense, this is not a fully self-consistent approach since the 
linearized equations are only appropriate when the turbulent velocities are small in the fluid 
frame, yet the dissipation of significant turbulent motions is the source of energy for the radiation. 
Nonetheless, we believe our approach should lead to a reasonable approximation to the truth. 
Simulations show that, when the only damping is numerical, the turbulence spectrum declines 
sharply toward shorter wavelengths. Thus, the short wavelength modes are legitimately in the 
linear regime, relative to the "equilibrium" background provided by larger amplitude, longer 
wavelength fluctuations, except that there exist non-linear couplings which cause the cascade 
of energy to smaller scales. A linear dispersion relation should at least provide a qualitative 
indication of the major effects. 



2.1. Photon Damping 

We first begin with a qualitative description of the different regimes of photon damping. 
When radiation pressure in a fluid is significant compared to gas pressure, momentum and energy 
can be transported by radiation in such a way as to damp out perturbations in the fluid. There 
are two relevant length scales: k^"^ = 1/neCrT (ng is the electron number density and ax is the 
Thomson scattering cross-section), the photon mean free path, and ~ k^^c/cg, the diffusion 
length (cs is the phase speed of long- wavelength acoustic perturbations). These two wavelengths 
define three characteristic regimes: 

1) Optically thin regime: When the wavenumber k > Inkx, photons can travel freely across a 
wavelength. The Doppler shift due to fluid motion creates a flux in the fluid rest frame that acts 
as a headwind for the electrons. As we will show later, this effect leads to a damping rate that is 
independent of k. 
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2) Non-diffusive regime: This is the range of wavenumbers < k < Inkx- In this regime, 
although a single wavelength is optically thick, photons can diffuse out of a fluctuation in a single 
wave period. This effect will prove especially important to compressive waves. 

3) Optically thick diffusive regime: When k < ki), photons are effectively dragged along with 
the fluid oscillations. Their diffusion can be described well by conventional transport coefficients 
( Weinberg 1972| ). If one thinks of the system as a single fluid, these correspond to shear viscosity 



and (a version) of heat conduction. Mihalas and Mihalas (1984) , and references therein, have 



derived these coefficients in the diffusion approximation. 

With these wavelength distinctions in mind, we now derive the exact dispersion relations for 
MHD modes damped by radiation transport. 



2.2. Radiation transfer equation 

I 1 

Mihalas and Mihalas (1984) derived the equations of radiation viscosity in the limit of 

time-steady, diffusive behavior. They additionally assumed pure absorptive opacity and LTE, in 
contrast to our assumption of pure isotropic scattering; however, this does not affect radiation 
viscosity. Our case involves time-dependent behavior and gradients that may be so sharp as to 
completely invalidate the diffusion approximation. Consequently, we must rederive the equations 
of radiation viscosity in a way that is appropriate for our regimes of interest. 

We write down the radiation transfer equation in a quasi-inertial "lab" frame which travels 
along with the local mean orbital velocity. We neglect rotation because we will be interested 
only in fluctuation wavelengths very short compared to a radius (in fact, for some purposes to 
make rotation negligible requires a stronger constraint to wavelengths very short compared to a 
disk thickness). We also neglect the thermal source function, absorption opacity, and stimulated 
scattering. The source function is then solely due to electron scattering. Evaluated in the lab 



frame and averaged over frequency, it is ( Pomraning 1973 , equation 6.1): 



1 /" ly f d(Tj^ 

5r(nf) = — / dvfdvid9,i^——{vi z^/,ni nf)/(ni,fj) (2) 

where /(n, u) is the specific intensity in the direction n at frequency ly, i and / subscripts indicate 
the initial and final photon respectively, gt is the Thomson cross section, and n is the direction of 
photon motion in the lab frame. If the fluid moves with velocity (3 (in units of c) relative to the 
lab frame, we have the following relations, correct to first order in /3: 



- = l-/3-(ni-nf), (3) 
^(i^i^z.y,ni^nf) = [l + /3.(ni-nf)]5[i.y(l-/3-nf)-z.,(l-/3-ni)]^, (4) 
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where the first relation is the familiar frequency shift due to Compton scattering, in which we 
have neglected terms of order /iiz/mgC^; the second is the transformation between frames for the 
Thomson scattering cross section, where we have made the approximations of isotropic scattering 
and negligible electron recoil. 

The first four moments of the frequency integrated specific intensity are: 

J = j;^Jdni{n) (5) 
H = -i- / dnnl{n) (6) 



4tt 

Kij = ^ J dQuiUjlin) (7) 

Lijk = ^ J dQuiUjUklin), (8) 

where dfi is s'mOdBdcj), n is the unit vector pointing in the {9,(j)) direction, and /(n) is the 
frequency integrated specific intensity. 

Integrating the source function (§) over solid angle and frequency and keeping only terms of 
order (3, we get: 

Srin) = (1 + 3/3 • n) J - 2/3 • H. (9) 

The full frequency-integrated radiation transfer equation in the lab frame, including only terms 
first order in /3 and neglecting emissivity and absorption is then 

1 ^^(") _^ . v/(n) = (1 + 3^- n)J- 2/3- H - (1 - n-;3)/(n), (10) 



ckx dt kx 

The last term on the RHS of this equation is due to electron scattering opacity, boosted from the 
fluid frame to the lab frame. This equation agrees with Psaltis and Lamb (1997), except that we 
have dropped terms second order in /3 and have ignored the temperature of the electrons. Taking 
the first moment of this equation (l/47r J dil.), we get: 

^ ^•^ + -iv-H = -/3-H. (11) 



c/ct dt kT 

Next, taking the second moment (l/47r J diln) gives: 

1 3H 1 

Finally, taking the third moment of the transfer equation, we find: 

Q^'' + -j^^^kLiju = — ■ H - Kij + PkLijk- (13) 

Now, to close these equations, we must make some assumption about the form of the radiation 
field. The Eddington approximation is equivalent to setting the quadrupole and higher moments 
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to zero. However, we want to consider the effect of radiation viscosity, which is only present 
if there is shear, and this requires a term of quadrupole order or higher in the radiation field. 
We therefore set all higher moments to zero, but retain the monopole, dipole, and quadrupole 
moments: ^ 

/(n) =Ii+l2n-nD + ^ [3(n • nqf - l] , (14) 

where ud is the direction of the dipole moment and ng is the direction of the quadrupole moment, 
and /i,2,4 are independent of n. Using this multipole expansion for the intensity, the fourth 
moment of the radiation field (equation |^ can be expressed in terms of the flux: 

Lijk = g [SijHk + SikHj + SjkHi] , (15) 

using the relation / {d^l/ATr)ninjnkni = {dij5ki + SikSji + 5u5jk)/l5. This result allows us to express 
the third moment of the transfer equation (^) as: 

i ^ + 5^ ^^''^ ■ " ^ + ^•'■^^^ = ^ ( J - 2/3 • H) - Kij + ^ [5ijP ■ H + PiHj + l3jHi\ . (16) 

Finally, we need to calculate the effect of the photons on the electrons. The rate of momentum 
transfer from the photons to the electrons via Compton scattering is: 

C = l dvdn'dm^0l^^n,^f{0){l - . n')^-^^^j^, (17) 

where primes denote the particles before scattering in the lab frame, is the electron velocity, /(/3) 
is the electron distribution function, and Ap is the momentum transfered during the scattering. 
We make the assumption that all electrons move with the fluid velocity /3, i.e. f{0) = 5^{0 - 0. 
In the limit of non-relativistic electron speeds, the momentum equation would be unchanged if we 
had instead averaged over a finite-width velocity distribution. We again assume the scattering is 
isotropic and we ignore terms of order hv/nieC? and higher. Performing the integral in equation 
(17) and keeping only terms of order /3, we find 

Ci = -^[PiJ + PjKji - Hi], (18) 

which is proportional to the RHS of equation (p^). We can then use this electron-photon 
momentum transfer rate in the fluid momentum equation. Ignoring all other forces, the fluid 
momentum equation is: 

where v = cf3. In most cases, C • /3 is negative, and thus there is generally a drag on the fluid due 
to collisions with photons. 



We assume that in the equilibrium state the radiation is uniform, time independent, and 
isotropic so that J = I. The higher order moments of the unperturbed radiation field are then 
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simply Hi = 0, Kij = 6ijJ/3, and Lijk = 0. We also assume that the unperturbed fluid is at rest in 
the lab frame. These assumptions greatly simplify the equations, and retain most of the physics 
of the waves. We assume that all perturbations vary with space-time dependence Q^O^-^-^t) ^ e.g. 
the perturbed mean intensity is J + 5 Je**-^ """""*^ ; for this reason we must also restrict attention to 
modes with fe/i S> 1. The perturbed radiation transfer equations are: 



lU! 



6J 



UJ 



6Hi 



— oKij 
krc 



ikj6Kij 4 

— -, + i:5l3iJ 

kx 3 

— J-[5»,k • (5H + ki6Hi + k,5Hi\ + ^5 J 
bkr ^ * J J »J 3 



5K, 



Solving these equations for 5 J in terms of 5(3 yields: 

4 



5 J 



-Jfl 



I UJ 



UJ (1 



I bj) H — k 
' 3 



1 



-I UJ 



k-5/3. 



(20) 
(21) 
(22) 

(23) 



where we have defined normalized variables k = ]<./kT, and uj= uj/kxc, so that the optically thin 
regime corresponds to |k| > 27r. The perturbed mean intensity disappears for modes with k _L 5/3 
since there is no compression of the radiation field. 



Next, the perturbed flux is: 



(5H = -J— 

3 (1 



(1 



I UJ) 



I UJ 



+ i fc^ 



5/3 - ik(k • (5/3) 



5 — 6i a; 



15tJ(l 



i ujY + i k^ (5 — 9i w). 



The perturbed collision integral is: 



5C 



-5 13 J - JH 



(24) 



(25) 



-r2 



For incompressive waves, in the limit that u:<^k <^ 1 (low frequency, optically thick limit), the 
momentum transfer rate is 

5C = --kTPradt 5(3. (26) 

5 



Thus, the photon-fluid friction is proportional to —k'^5(3, which looks like a S/'^v term, with a 
constant of proportionality ^PradkT /5c. This is the same as the radiation viscosity term derived by 
multiple authors, e.g. Mihalas and Mihalas (1984). Indeed, the photon viscosity computed by Loeb 
and Laor (1992) for a steady shear flow gives exactly this viscosity term, except with a factor of 
10/9 from considering the exact differential Thomson cross section (rather than assuming isotropic 
scattering as we did above). Including polarization introduces another small correction factor ( Hu 



and Sugiyama 1996 ). In the large k/kx (optically thin) limit, the friction approaches a constant, 
which is what one would expect for electrons in a uniform radiation field; since the wavelength is 
much shorter than a photon mean free path, each electron sees the averaged radiation from several 
wavelengths. The friction changes when a; 7^ to take into account time-dependent diffusion of 
the radiation as the wave oscillates. 
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2.3. MHD equations 

With the perturbed radiation quantities in hand, we now write down the perturbed MHD 
equations of motion. We define the z axis as the direction of the magnetic field, and ignore both 
gravitational potential gradients and rotation, in keeping with our restriction to modes with 
kh S> 1. Ignoring rotation is valid for ^ 17, or fc/i S> Cg/vA- By omitting rotation effects, we 
restrict our attention to wavenumbers short enough that the Balbus-Hawley instability does not 
operate. The effects of vertical gravity on radiation waves in accretion disks has been considered 
by pammie (1998) , who found unstable "photon bubble" modes. Our equations ignore these 



modes; however, we do include radiation viscosity, which Gammie (1998) ignored. The MHD 
equations are 




1 



(B • V)B - C = 



'dt 

where Pgas = NksT and /c^ is Boltzmann's constant. 



dp 
dt 

- V X (v X B 



+ V • (pv) = 




(27) 



Assuming an equilibrium state with p and B constant, and v = 0, the perturbed equations 



are 



-iujpSv + iMPgas + -r [k(B • SB) - (B • k)6B] - 6C = 

6p k • 

p UJ 

u)5B = B(k • 5v) - (B • k)Jv 



(28) 



We will assume that 5Pgas/5p = Cg is constant in what follows. Combining these equations gives: 



(k-<5v)k- , 

c y \ c J 



2 ^ 



(k • Jv)k+ k1 6\r— kz (k • 5v)z— kz 5vz\s~ 



pkxc 



5C = 0, (29) 



where va = \/i?^/47r/?, the Alfven speed. 



Now, define k = /c(sin0x + cos 0z). Setting the determinant of equations (29) to zero, and 
using equations (p3) and (p5|), we find the following dispersion relation: 



Al I (Al + ^3 cos ef + {A1A2 - Al) k^ 
where we have defined the auxiliary quantities: 

{up' — v\^^ cos^ 9 + i tUr") Di — i cUF {I — iuJ) 







(30) 



^3 



-(Cg + v'a)Di- ujT (5 - 6i uJ)(l - i aJ)/L>2 
cos ODi 
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Di = {l-iLjf + - 

5 

L>2 = Ihuj {l-iujf + i 1? (5 - 9i TU) 
VA = va/c 

r = T/kTC = ^Prad/pC^ = ^cl/c^ (31) 

In the limit Prad = 0, this dispersion relation becomes the MHD dispersion relation for Alfven 
modes and magnetosonic modes (e.g. [Jackson 1975| ). For nonzero Pradj this dispersion relation 
admits a variety of modes: modified versions of Alfven modes (incompressive) ; fast and slow 
magnetosonic modes (compressive); and radiative (electromagnetic) modes. On account of the 
complexity of the full dispersion relation, we will discuss some simplified limits. The dispersion 
relation for Alfven-like modes factors out separately, Ai = 0. We will discuss this branch in the 
next subsection. When = Oor0 = 7r/2, the part of the dispersion relation in brackets simplifies 
considerably; we will look at the compressive modes for these propagation directions in § ^.5| . 



2.4. Dispersion relation for incompressible waves 

The case of = yields the modified Alfven modes, for which k • (5v = and 5vz = 0. 
Since these modes are incompressible, 5J = Q and (5H simplifies drastically. The equation = 
becomes: 

(32) 
(33) 



up' — k^v\ cos^ 9 + iAiU) = 0. 



where 



A, = r- 



-i uj — uj'^ 



{i-iujy + l k 

If VA ^ c, then | ZJ \ <^k (since ~ kvA) and we can expand: 



A. 



5+ k 



5i io 



5- k^ 

5+ k^ 



(34) 



Substituting this expression into equation (|32|), we can solve for uj: 



'-\ 



k\\ cos2 OAk- r 



.5+ k 



2 -2 
I .2 .-T fe 



■2 ' 



(35) 



where Ak = {1 + 5 T (5— k )/(5+ k )'^) ^. This dispersion relation is valid for any T and for 
VA ^ c. For T<^VA, it simply describes damped Alfven waves, and it agrees with pedamzik et al 



(1998) in the diffusion, optically thin, and overdamped limits. Note that the damping rate we 



have computed bridges the optically thin and thick limits, which have been examined separately 
by other authors (Subramanian and Barrow 1997, Jedamzik et al. 1998). 
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We also find higher frequency modes with \uj\ ^ T, illustrated in figure ||. These modes are 
non-propagating in the optically-thick limit and damped on a timescale ~ (/c^c)"^, but travel 
at c/\/5 when their wavelengths are shorter than a photon scattering length. Their dispersion 
relation is to first order the solution of Di = 0, which means that 5f3 — > 0, so that these modes 
are simply electromagnetic. The speed of these modes is likely an artifact of the particular form 
taken for the moment hierarchy closure; their speed increases when we include moments higher 
than quadrupole. 



2.5. Dispersion relation for compressible waves 

Next, we consider the opposite limit of the dispersion relation: strongly compressible waves for 
which k||(5/3. This condition permits two sorts of waves: ^ = (k||B) which is simply a radiation 
damped sound wave; and 9 = 7r/2 (k _L B), the fast magnetosonic wave. The slow magnetosonic 
wave disappears when we impose k.\\5f3 since it has a velocity component perpendicular to the 
wave vector. The damping rate of the slow magnetosonic wave is intermediate between the Alfven 
wave and the fast magnetosonic wave, which is why we do not treat it here. The fast wave also 
has a small k _L (5v component when it propagates at an angle < 6 < 7r/2; we ignore these waves 
here since they will have damping rates in between the 6 = and 6 = ir/2 cases. The dispersion 
relation for the 6 = 7r/2 case is: 



UJ 



(36) 



where we have defined 



A, 



Z UJ 



^(5-6i ..) 



In the diffusive limit, k <^ ki), the fast magnetosonic dispersion relation becomes: 



1 + 3% 



k'^icl + v\ + 4) +iTui 



2 -2 k^ 2' 
5 9 uj"^ 



+ 



1 

kT 



0. 



(37) 



(38) 



where = APrad/'^P is the sound speed due to radiation, and we have assumed that all velocities 
are smaller than c. Since k <^ k£), k is very much smaller than kx, so it makes sense to expand in 
terms of k/kx- In terms of this expansion, the zeroth order solution to equation ( |38[ ) is 



UJo 



vi + ci 



\ 1 + 



3c2 



k. 



(39) 



We now substitute this approximate solution into the next higher order term in equation (I 
obtain a quadratic equation for uj with the solution: 



and 



UJ 



±U)Q 



QkT{l + R?{l + -^R) 



W + -{l + R) + Q{- 



) n R + 



5 ' 



(40) 
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Fig. 1. — Dispersion relation for k _L 5/3 waves with va = Cr = 0.3c and cosO = 1. The sohd hnes 
are Alfven modes, while the dotted lines are the electromagnetic modes. The dashed line is the 
damping rate from equation (35). 
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2 2 

where R = 1/ T= — = ^-nd we expect i? ^ 1 in all applications of interest here. In the 



limit Cg = = 0, this equation has the same form as equation (52) in Peebles and Yu (1970) . In 
the limit that all velocities are smaller than c, this equation becomes: 



ck^c 



2„2 



Thus, in the diffusive limit, the damping rate has the usual oc k"^ dependence. 



In the optically thin limit {k ^ kx)-, the dispersion relation becomes: 



u; = ±k^cl+v\-"f, (42) 

giving the same damping rate as for incompressible waves in the optically thin limit. Compressible 
waves are damped at the optically thin rate when the time for photons to diffuse across a 
wavelength is comparable to or shorter than the wave oscillation period. We do not have an 
analytic formula for the dispersion relation for ko ^ k <^ kx, so the full dispersion relation must 
be solved numerically, as shown in figure ^. A more accurate formula for kjj can be found by 
equating the diffusive and non-diffusive damping rates: 



kD = y'3(3c2 + c2 + v\)kT/c, (43) 

which is where the compressive damping rate reaches its near-maximum. 

When pc^ S> Prad S> Pgas (i-s. Cr <^ c, Cg <^ Cr) and VA *C Cr the damping rate for 
magnetosonic waves at small k is ck'^/dkx- The damping rate for Alfven modes, however, for small 
/c is ~ r A; /lO = 3{cr/c)'^ck'^/{10kT)', i.e. it is smaller by a factor ~ (crjcf' . Thus, magnetosonic 
modes are damped much more strongly than Alfven modes in the small A;, small F limit. The 
reason is that compressional waves continually compress the radiation field, which diffuses out of 
the wave, causing the wave to lose its pressure support, and thus damping it out. Since there 
is no compression in the Alfven modes, photons only diffuse out of the wave perpendicular to 
5v, creating a quadrupole moment in the radiation field which leads to viscous damping, and 
that is much weaker than diffusive damping. A comparison of the dispersion relation for Alfven 
modes and magnetosonic modes is shown in figure |2| with va = = 0.1c, = (r= 0.03). In 
the optically thin limit, both waves are damped at a rate r/2, since the radiation is isotropic 
and uniform and thus the damping is just due to the dipole moment of Doppler shifted photons 
in the fluid frame. The phase speed of magnetosonic waves with these parameters is somewhat 
greater in the diffusion {k < k^) limit than at larger k because the radiation adds to the restoring 
force. As a corollary, the waves are mildly dispersive for k ^ ko- Also shown are the analytic 



approximations to the damping rates, equations (p5|) and (40) 



More complicated behavior can occur when va *C Cr because the removal of radiation pressure 
support when k exceeds k^ is a relatively more important effect. In figure |3|, we show the 
compressive dispersion relation (equation ^) for va = 0.001c, Cr = 0.01c, and Cg = 0. There are 
three separate cases to consider for the fast magnetosonic modes: 



- 13 - 




Fig. 2. — Plot of the real and imaginary parts for Alfven (solid lines) and magnetosonic (dotted 
lines) with va = Cr = 0.1c, Cg = 0. The dot-dash line in the top panel is equation (|^); the dashed 
line (which overlaps the solid line) is equation (35). 
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1) In the diffusion limit, the modes are simply sound waves supported by radiation pressure 
since the diffusion timescale is much greater than the wave period. 

2) In that portion of the non-diffusive regime in which ko < k < T/{2va), the drag due 
to radiation escaping from the waves is greater than the magnetic restoring force, so the wave 
becomes overdamped. This is the region with Re{uj) = in figure ^. 

3) For k > V/[2va)-, the radiation isotropizes on a timescale much shorter than the wave 
period. In this range of wavenumbers, fast magnetosonic waves propagate with phase speed 



^ + c^, but damp due to radiation drag. 

There are additional overdamped modes which occur when the velocity perturbation is so 
small that the fluid reaches a terminal velocity and thus the mode damps out before it can 
oscillate (Subramanian and Barrow 1997). In the diffusive limit when Re{co) = 0, the radiation 
always has time to isotropize, so the photon damping rate is the same as in the optically thin case. 



The magnetic restoring term in equation (29) with kz = and k||5/3 balances the optically thin 
radiation collision term when 

uj = -I— „ ^ ^ . (44) 

This agrees exactly with the damping rate for the lower dotted curve in the upper panel of figure |^. 
For k ^ O.lkT, the upper dotted curve is the solution of D2 ~ (for very large A;, w ~ —i5kTc/9), 
which means 5/3^0, so that this mode becomes electromagnetic (but non-propagating) in the 
optically thin limit. 

We also find propagating electromagnetic modes (not plotted) whose velocities approach 
and c-\/3/5 in the optically thin limit. These speeds are again artifacts of our closure relation. 



2.6. The nature of the transport 

Why is it that photon transport is so much more effective in damping compressive waves than 

incompressible ones? One way to understand the contrast is to take a closer look at equations ( p4[ ) 

and (^). In the incompressible case (k • 5/3 = 0), (5H ~ {A/3)JS(3 when uJ and k are both <C 1. 

That is, the perturbed flux is simply 4/3 times the mean intensity shifted by /3. However, the 

force felt by the electrons is the difference between 6H and (4/3)J(5/3, so the two very nearly 

cancel. The remainder after this near-cancellation is the retarding force due to the photon shear 

2 

viscosity, and has magnitude ~ (uJ -|- A; ) J. However, the relative importance of the photon drag 
is characterized by the ratio F, the ratio of the photon inertia to the fluid inertia, and this is <C 1. 

On the other hand, in the compressible case, there is an additional contribution to 5H 
(and therefore 6C) that is, in the very low frequency limit, ~ (4/3) J5/3, much greater than the 
contribution of shear viscosity when ZJ and k are <C 1. This new contribution is the result of 
photon diffusion. In a single-fluid picture this can be thought of as a sort of thermal conductivity 



-15- 




Fig. 3. — Plot of the real and imaginary parts for magnctosonic waves with va = 0.001c, Cr = 0.01c, 
Cg = 0. The dotted lines are overdamped waves, while the solid lines are magnetosonic waves. 
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( Weinberg 1972| ), but with respect to a peculiar equation of state (cf. equation 23). It damps the 
waves much faster than shear viscosity because the diffusing particles in this case are exactly those 
responsible for the wave's restoring force. 

In principle, photon-electron scattering could also lead to a magnetic diffusivity by 
creating a new source of electrical resistivity. In the optically thin limit, the resistivity is 
rjr — 1.4 X 10~^^ (cr/c)^ s. In most parameter regimes this photon resistivity is small compared to 
the resistivity due to ordinary electron-ion Coulomb scattering, r/p ~ 1.4 x 10~'^T~^^'^ s; however, 
in the corona, where Tc ~ 10^ K and Cr ~ 0.1c, photon resistivity may compete with Coulomb 
resistivity. Due to flux-freezing, the damping of the turbulent motions also damps the magnetic 
field fluctuations; however, due to the small resistivity, magnetic flux is still conserved. 



3. Applications to conventional accretion disks 
3.1. Context 

In the preceding section we characterized the effects of photon diffusion and viscosity in terms 
of the rate at which they cause damping of linear MHD waves. In this section we will evaluate 
how effective these processes may be in dissipating fluctuations in accretion disks. To gain a sense 
of scale, we begin this section by estimating the corresponding damping rate for several other 
proposed dissipation mechanisms. Although the natural unit of time for the dispersion relation 
was the photon scattering time, in the context of disks the natural unit is (the inverse of) the 
orbital frequency fi, so we will quote all rates in that unit. As a further set of reference rates, 
in this subsection we will also establish the relevant standards of comparison for several different 
questions of interest. 

Ordinary molecular viscosity (due to ion-ion collisions) creates a damping rate 

T^oi = Ukhf^^n, (45) 

-J (^coU Cs 

where h is the (half) disk thickness and acoii is the collision cross section. If acou is the Coulomb 
cross section, for example, (Jt/cTcoU ~ {ksT / meC^)"^ / In A, where In A is the usual Coulomb 
logarithm, ~ 30. Tmoi/^ is generally a very small number. Ordinary viscosity is rendered even 
less effective because the magnetic field suppresses transport perpendicular to the field. 

Transit-time damping (and the associated Landau damping) has been suggested by Quataert 
(1998) as the dissipational mechanism in accretion disks, particularly when the ion temperature 
is much greater than the electron temperature as in advection dominated accretion flows. As a 
fiducial point, we quote its rate (as calculated by Quataert 1998) for a single-temperature plasma: 

2 2 

r,,, ^ 0.2 cos e sini 9 {kh)l ^ ' (^) ' n, (46) 
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where is the ion Larmor frequency. When k is greater than an inverse ion Larmor radius, 
oc {khy/'^. Unless the magnetic field is exceedingly strong, ^l/Qi ^ 1, so that r^^^^ ^ fi. 

Depending on the question being asked, any candidate damping rate should be compared 
to one of three fiducial rates: the growth rate (absent dissipation) of the MHD waves (as, for 
example, due to magneto-rotational instability as in Balbus and Hawley 1991); the inverse time 



for waves to cross a disk scale- height; and the "nonlinear frequency" or inverse "eddy turnover 
time," the rate at which energy moves between modes due to nonlinear coupling. 

If the damping rate exceeds the non-dissipative growth rate, the fluctuations are unable to 
grow at all. This is a strong statement, for the growth rate of the magneto-rotational instability is 
generally ~ 17. 

When the damping time is short compared to a disk scale-height crossing time, waves cannot 
carry significant energy from the midplane to the disk surface (see §4.2). This, too, may require 
very rapid damping, for the wave crossing time can be as short as ~ (for diffusive regime fast 
magnetosonic modes). The time for other modes to traverse a disk thickness is somewhat slower: 
~ {Qvyi^/cs)"^ for pure Alfven modes, ~ {Q^v\ + c^/cs)'^ for fast magnetosonic modes with 
k > ko and uj ^ 0. 

Thirdly, as emphasized by Gruzinov (1998), exceeding the nonlinear frequency at some 
wavenumber is the relevant criterion for deciding whether the damping can cut off the "inertial 
range" of turbulence at short wavelengths. The nonlinear frequency is defined by uJni{k) = e/kEk 
where Ek is the energy density per unit wavenumber in the turbulent spectrum. The rate of 
energy dissipation per unit volume, e, is determined solely by the accretion rate (equation 1) and 
disk height, while the total energy in fluctuations may be related to the accretion rate if we know 
the ratio between the trace of the fluctuations' stress tensor and its r — component (under the 
assumption that it is this last quantity which accounts for angular momentum transport in the 
disk). That is, the total energy density in fluctuations, volume averaged, is 

1 i-kmax Tr(T) MVt 

\Tr{T)= dkE, = i^^, (47) 



where Tij = {p{5vi5vj — 5vA,iSvA.j)) , and () denotes volume averaging (iBalbus et al. 1994| ) 



Suppose the fluctuation spectrum is a power-law E^ oc k " from kmin = ir/h to kmax- Then, 

n. (48) 



Trip 



{n - 1) Tr{T) 
If kmax > kmin and n> 1, then 



- 



vr y \k 



kh\ 



^max 



n-1 



In the simulations of Brandenburg et al. (1995) and Stone et al. (1996), Tr,j,/Tr{T) ~ 0.1, very 
roughly. 
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3.2. Radiation-pressure dominated disk 

First consider radiation-pressure dominated disks. In this case, the Shakura-Sunyaev solution 
(in which Trij, is set equal to ap) yields two important results about the equilibrium. The disk 
aspect ratio is 

^ = | = , (50) 

r 2 X 

where m is the accretion rate in Eddington units (for unit efficiency) and x = rcP/GM is the 
radius in gravitational units. We have ignored all relativistic factors. In addition, the (half) 
optical depth is 

2c 4x1 
mm Sam 

where, in consonance with the result of simulations, we have ignored any contribution of magnetic 
pressure to disk vertical support. 

We may now use these facts to evaluate the rate of photon damping. In this case, 

— = -mx 2 . (52) 

Typically c^ <C c (r<C 1) in thin accretion disks. For k < ko and radiation pressure larger than 
magnetic or gas pressures, compressive modes damp at a rate 

r-p = ^^kh)'n. (53) 

Incompressible modes damp more slowly at small k: 

"pmc ™ 2 

±d ^[[^ C54N 

So only for large m and very small radii will the incompressive damping rate equal or exceed the 
compressive. In the optically thin {k > lirkr) limit, both damping rates are constant and equal to 

Ythin ^ 3q,-10. 

Thus, we immediately see that compressible modes damp extremely rapidly. At the longest 
wavelengths the damping time is, not surprisingly, the same as the thermal time, {aQ)~^. They 

are, in fact, the same process — photon diffusion out of a region ~ /i in size. These modes damp 
so quickly because in this regime photons provide most of the pressure; consequently, it is their 
diffusion rate, not the ions', which controls the damping rate. 

To sec just how rapid the photon damping is, we may compare it to, for example, the rate of 
transit-time damping. In this context of radiation-dominated accretion disks, 

Tttd - 6.3 X IQ-^^ cosOsin'^/^ e{khf/^{cg/csf/^M~^^^n, (55) 

where Mg is the mass of the central black hole in units of IO^Mq. Because Cg/cg ^ 1 when 
radiation pressure is dominant, Fud is very slow indeed compared to even F^"*^. 
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So long as a < 1 , the damping rate for compressible modes does not exceed the Balbus-Hawley 
growth rate. However, the damping rate may well exceed the nonlinear frequency even for the 
longest wavelength modes, for 



-pcomp 
tA 



36 



i)rr(r) (kh\ 

vr / 



P 



(56) 



This expression follows from the fact that the Shakura-Sunyaev parameter a = 
\Tr(j>/Tr{T)]\Tr{T) /p]^ where p is the total pressure. If the spectrum has the Kolmogorov 
slope (n = 5/3), the photon damping rate is greater than the nonlinear frequency for wavenumbers 
not much greater than vr/Zi unless Tr(T)/p very much less than one (in the simulations of Stone 
et al. 1996 this quantity was ~ 0.01 - 0.1). 

Even if photon damping does not overcome the fluctuations at longer wavelengths, it is still 
likely to end the inertial range of turbulence. The maximum photon damping rate (the optically 



thin limit) is achieved at k^h 
rate to UJni, we find 



Ycomp 

d,max 



6/q, where the damping rate is F 



comp 
d.max 



3^}/a. Comparing this 



in 







71 — 3 


■Tr{T)' 




Tr{T)_ 




. P . 



n n-2 



(57) 



So long as n < 2, it is almost guaranteed that > ujni at some wavenumber. For instance, for 



n 



5/3, we find that 



-pcomp 

d,max 



43 



10- 



loorr(r) 



'Tr{T)\ L P J ' ^^^^ 
where we have normalized to fiducial values in the ballpark of what is seen in simulations. We also 
emphasize that this equation is independent of x, m, and mass of the black hole. Thus, radiation 
damping of compressive modes can be quite strong whenever radiation pressure dominates. 

If there is any significant azimuthal field, i.e. va4) — Cs, all long wavelength (i.e. k not too 
much larger than 0,/va) modes are at least partly compressible (Blaes & Balbus 1996). Simulations 
indicate that va ^ Cs, so the Balbus-Hawley unstable modes are very nearly incompressible. 
However, there can still be significant coupling between incompressible and compressible modes. 
In simulations of nonlinear magnetohydrodynamic turbulence by ptone et al. (1996) , the density 
fluctuations {6p'^)^/'^/p ~ 5 — 8%, while the velocity fluctuations {6v'^)^/'^ / Cg ^ 15%. Since 
6p/p = k • 5v/uj ~ k • 6v/cs, then k • (5v ~ 0.5. Thus, the waves in the turbulence spectrum have 
a rather large compressive component. Also, in these simulations, pressure waves are seen which 
are not present in the non-turbulent state, indicating that the turbulence does create compressive 
waves (John Hawley, private communication) . Another way to quantify the fraction of compressive 
turbulence is to take the power spectrum of the vortical and compressive components of the 
velocity, v = Vvort + Vcomp such that V • Vvort = and V x Vcomp = 0. MHD shearing box 
simulations by Brandenburg et al (1995) show that the power spectrum amplitude of Vcomp is 
about 10% of that of Vyort ( prandenburg 1998| ). 



- 20 - 



In addition to the canonical solution for radiation-dominated disks, there are also several 
extensions of this solution to which photon damping is relevant. At high accretion rates, rh > 1, 
the radiation pressure causes the disk to puff up, creating a "slim disk" ( [Abramowicz et al. 1988 ). 



Slim accretion disks have larger luminosities than thin accretion disks, making the effects of 
radiation damping much stronger. The standard thin disk equations cannot be applied to slim 
disks, since some of the radiation is carried radially inward through the disk rather than being 
radiated locally. However, the slim disk solutions look similar to the Shakura and Sunyaev (1973) 



solution in the limit of large rh ( Szuszkiewicz et al. 1996| ), so we expect our criterion for radiation 



dissipation to apply, under the assumption that slim disks can be approximately described by the 
thin disk equations. 

Radiation pressure-dominated disks in which T,.^ = aP^ot are viscously and thermally unstable 
( Shakura and Sunyaev 1976| , Lightman and Eardley 1974| ). To cure this, some have suggested that 
the viscosity is proportional to the gas pressure rather than the total pressure, i.e. Tr^ = aPgas- 
Indeed, as we will discuss in §5, photon diffusion may decouple the radiation pressure from the 
MHD fluctuations, leading to just this sort of result. If so, such disks effectively have a much 
smaller a, and consequently radiation viscosity is much more efficient at damping perturbations 
since 1) the turbulent velocities are much smaller, reducing the nonlinear frequency and 2) the 
optical depth of the disk is much larger, so the radiation pressure at disk center is larger. Since 
the disk is still supported by radiation pressure, the height is the same as for aPtot disks, and thus 
r also remains the same. However, the non-linear frequency is reduced by a factor of Pg/Ptot, 
which is given by 

^' - (59) 



Ptot 27 aTiJi 

In these disks, even the incompressible damping rate can beat the nonlinear frequency at A; ~ kx, 
the wavenumber at which T^^^ reaches its maximum value. For example, for n = 5/3 and 

Tr{T)/Ptot = 10, 

— ^— - = ^Tifa^x-'^. (60) 

The optical depth in these disks is given by 

r = 5 X lO^a'tmtMg^x"!. (61) 

Using this expression for r, we see that the incompressible damping rate beats the nonlinear 
frequency out to a radius of 

14 ^7 16 , X 

X < 178mT9Mg5^a"57. (62) 

Since the compressive damping rate is always greater than the incompressible, radiation damping 
will be important for a large range of radii if the viscous stress scales with gas pressure rather 
than radiation pressure. 



Our discussion of accretion disks so far has neglected the vertical stratification of density and 
radiation pressure since we have been using values computed from a one-zone model. Because the 
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radiation damping rate is proportional to the radiation pressure, we would expect the damping 
to be relatively more important in the interior of the disk, so that more dissipation occurs deep 
inside the disk (provided the nonlinear frequency is independent of height). 



3.3. Gas pressure-dominated disks 

The damping criterion we have discussed only applies to the P^ad 3> Pgas regions of the 
accretion disk. When Pgas or Pmag ^ Prad, the damping rate becomes (cA;^/2A;t)c^/(3c^ + + 
(cf equation pO| ) in the diffusive regime, so radiation damping will not compete with the 
nonlinear frequency. Thus, we expect that the radiation damping will only be important in the 
radiation-pressure dominated part of an accretion disk. The radius at which radiation pressure 
equals gas pressure is given by: 

Xtrans = ISSiaMs)^ m^^ {I - /)^ (63) 

where Mg is the black hole mass in terms of W^Mq, f is the fraction of energy lost to a corona, 
and we have assumed Thomson scattering opacity. Since this radius is very insensitive to the black 
hole mass, we expect radiation dissipation to be important in the range of radii in which most 
of the luminosity is created for black hole X-ray binaries, Seyfert galaxies, and quasars. There 
is a rather strong dependence on the luminosity relative to Eddington, so radiation dissipation 
won't play a role for objects with small rh. We have computed a disk model which includes both 
radiation and gas pressure, and compared the nonlinear frequency with the numerical root of the 
dispersion relation for compressive waves. We find that the radius at which the damping rate 
exceeds the nonlinear frequency (for some k) is typically at Pgas — a few xP,.ad- 



3.4. Growth of the radiation field 

So far our discussion of accretion disks has assumed that they are already radiating. Since the 
radiation is derived from the dissipation of kinetic energy into electron thermal energy or photon 
energy density, we have only showed that radiation damping provides a dissipation mechanism 
which gives a self-consistent disk solution. Another stronger question to ask is: if a disk is in 
a state in which radiation pressure is small relative to gas pressure, for what parameters will 
the radiation dissipation cause growth of the radiation field, causing the disk to find a radiation 
pressure-dominated equilibrium? 

The rate of change of the radiation field is given approximately by: 



dPrad Q 

— - — = —mm 
dt 2h 



1,(^ 



max 



cPrad 



h{l+T) 



where the first term on the right hand side is the rate of creation of radiation due to photon 
dissipation of turbulence {Q is given by equation 1), and the second term is the rate of escape of 
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radiation from the disk. Now, Q ~ cPrad/{^ + 't) in equilibrium, so a steady-state radiation field 
can only be achieved for {Td/tjOni) ;^ 1. In general, the maximum damping for compressive waves 
r/2 occurs for k ^ kf) kxCs/c. For a general disk, the maximum ratio of radiation damping to 
nonlinear frequency is given by: 

(65) 

where Vg = GM/c^, and the constant of proportionality is of order unity. For incompressive 
modes, the maximum damping rate occurs for ~ A:^, where 

i^ir^] rhx-lr'--^. (66) 

These expressions are valid for optically thick disks which may be radiation pressure or gas 
pressure supported, and have Q given by equation (1). Whether either of these is greater than 
unity depends on what state the disk begins in. We consider one such starting state in the next 
section: an advection-dominated disk. 



4. Unconventional accretion disks 

4.1. Advection-dominated disks 

In an advection-dominated disk, the equilibrium depends on the fact that the cooling 
timescale is much longer than the accretion timescale and thus the heat is advected inwards rather 
than being radiated locally. If radiation damping is strong enough to cause growth of the radiation 
field, then the radiation will damp out the turbulence and most of the heat will go into radiation 
rather than proton thermal energy which gets advected. 

To estimate when radiation pressure is subject to growth, we assume a steady state disk with 
electrons of a constant temperature in which the viscous stress is generated by magnetic fields 
which create a turbulent cascade to smaller wavelengths. Using the criterion of Narayan and Yi 
(1995) for the existence of an advection-dominated solution (m < 0.5a^), we find r < (a/0.1)x~^/^, 
so advection-dominated disks are usually in the optically thin regime. Since the disk is optically 
thin, the radiation damping is given by F^^*" for either compressive or incompressive modes. 
Comparing the damping rate to the slowest nonlinear frequency (at kh = tt), we find that the 
criterion for radiation growth using equation (||) is x < 0A5[a^Tr{T)/Tr^]'^/^ , which means that 
radiation viscosity will not cause optically thin advection dominated disks to cool and radiate. 
This also means that whatever radiation is produced in an advection-dominated accretion flow 
cannot be produced by the photon damping mechanism. 
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4.2. Corona-dominated accretion disks 

An alternative disk equilibrium has been proposed by Svensson and Zdziarski (1994), in whicii 



all of the angular momentum transport occurs within the accretion disk, while a fraction / of the 
associated heat released occurs above the disk in a corona. Their equilibrium relies on the idea 
that the energy can be efficiently transported from the disk to the corona somehow, presumably 
through magnetic or acoustic waves. Since there is no outgoing radiation flux within the disk, 
its equilibrium density and gas pressure are much greater than in the radiation-supported case. 
However, unless / is very close to unity, there will still be a signiflcant region of the disk in which 
radiation pressure dominates (see figure 2 of Svensson and Zdziarski 1994). In this case equations 
( |53| ) and ( |5^ still apply, and the radiation damping time for compressive waves is less than the 
wave crossing time, 2TrQ^^ for 

For incompressible modes, the crossing time 2ttCs/{0,va) is greater than the damping time for 



Sx^VA 
rn?a{l - f)cs 



^^;^\/-XT^rAT:-- (68) 



Thus, only a limited range of wavelengths can successfully carry energy to the corona. 

When gas pressure dominates, if va — Cg, the crossing time for hydromagnetic waves is about 
2tt/^}. For incompressive waves with wavelengths less than the mean free path of a photon, or for 
compressive waves with k > kf), the damping rate will be F^''*"'. For / — 1 disks, 

knh = 2 X 10^mix"ia"^(l - //2)"im|, (69) 

while 

kTh = T = 1.8 X 10Sm3/4x-3/4a-iM8/^(l - f /2)-^/^{k/ kt)- (70) 

Now for waves to be damped by photon viscosity before they can escape from the disk requires 
T/2Q > 1. This ratio is: 

neglecting relativistic factors, where T and p are the density and temperature inside the disk, 
and we have used the equations from the appendix of ^incell and Krolik (1997") to evaluate the 



disk parameters for / — 1. Thus, only extremely short wavelength waves may be damped rapidly 
enough, and then only in rather extreme conditions (relatively large m and small x). If VA/cg <^ 1, 
the requirements for damping incompressible Alfven waves may be relaxed somewhat, but unless 
this ratio is very small, the qualitative conclusion is unlikely to be altered. 
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5. Discussion 

We have shown in the previous sections that the effectiveness of radiation in damping fiuid 
motions depends strongly on the ratio of radiation to gas pressure. As noted in equation (]63|), 
radiation tends to be most important in the inner parts of accretion disks, which are, of course, 
the most important for energy release. At least some part of the disk is radiation-dominated when 

m > 1.0 X 10-3^^iA6^-i/8^^-i/8^ (72) 

where Xmin is the inner radius of the disk. If the central object is a black hole or a weakly- 
magnetized neutron star, we may expect Xmin to be the radius of the marginally stable orbit, = 6 
in the limit of a spinless black hole, and ^ 1 as the spin of the black hole approaches its maximum 
possible value. However, if the disk does not extend in so far, whether because the central mass 
is a strongly-magnetized neutron star, or a larger object such as a white dwarf, the minimum 
accretion rate for which at least part of the disk is radiation-dominated rises, and may become 
impossibly high. 

The remainder of this section, in which we outline the consequences of radiation damping in 
accretion disks, is divided according to consequences applicable to radiation-dominated disks and 
those applicable to the gas pressure-dominated case. Whether one set or the other is relevant to a 
given disk depends on how it fares according to the criterion of equation (|7^). 



5.1. Radiation pressure-dominated disks 

Two qualitative physical consequences follow from the strength of radiation damping in 
photon pressure-dominated disks. First, dissipative heating is delivered to the electrons and 
photons through radiation scattering, and not to the ions. Because it is the electrons that cool the 
gas through the creation and upscattering of photons, the only energy exchange process involving 
the ions is Coulomb scattering. This mechanism should keep the ion temperature very close to the 
electron temperature. If the average energy of photons is less than /3^meC^/3 + 4A;BTe, where is 
the electron temperature, then the photons will receive most of the energy from scattering (Psaltis 



and Lamb 1997). The P'^rUeC^ term represents a modification of the Compton temperature due to 



bulk Comptonization. 

Second, the process by which these disks shine may be thought of as a sort of "bootstrap": if 
the disk were initially free of radiation, any initial photon creation by the electrons would lead to 
wave dissipation that heats the electrons, and therefore leads to more radiation. The question of 
what makes near-Eddington accretion disks shine has a tautological answer: bright accretion disks 
shine because they are so bright. 



That MHD fluctuations should be present at all is likely due to the operation of the 
magneto-rotational instability identified by Balbus & Hawley (1991). This instability grows at a 
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rate ~ kvA for wavenumbers k < \/^^/va when the magnetic field is weak (i.e. va < Cg). The 
compressibHity of the growing modes is shght, so the corresponding radiation damping rate should 
be a fraction of the pure compressive rate, as given by equation ([5^). If most of the torque in the 
disk is due to magnetic fluctuations, the ratio between the magneto-rotational growth rate and the 
radiation damping rate is then at least ~ 10{cs/vA)/ikh). We therefore expect the linear growth 
of MHD fluctuations to proceed unaffected by radiation damping. 

However, shorter wavelength waves are not amplified by the magneto-rotational instability. 
Instead, they are pumped by nonlinear coupling with the longer wavelength, growing modes. 
Because the radiation damping rate is oc k"^ in the diffusive regime, compressive modes excited 
by nonlinear coupling will be strongly damped. In other words, provided only that the nonlinear 
coupling between incompressible and compressible modes is reasonably strong, the "inner scale" 
of the MHD turbulence will be not much shorter than its "outer scale." Any turbulent "inertial 
range" will be severely limited. 

This fact leads to several other results. At a purely technical level, if short wavelengths are 
all severely damped, the life of the numerical simulator is made much easier, for there is no need 
to strive for very fine spatial resolution. 

More physically, radiation damping may play an important role in regulating the value of the 
"viscosity" parameter a. The magnetic part of the stress causing angular momentum transport 
may be written in the form 

Tr<p = ^J d^k5Br{k)SB;{k). (73) 

If there is little power in the fluctuations at wavenumbers much more than ~ the angular 
momentum transport is reduced below what it would otherwise be. Disks in this situation 
would then maintain rather larger surface densities. Increased optical depth also leads to greater 
radiation pressure for fixed emergent flux. 

Another consequence for turbulence in radiative disks is that the ratio of the sound speed 
to the Alfven speed changes with wavelength. In the diffusive regime, Cs ^ Cr ^ va, leading 
to a large plasma (3 = Ptot/Pmag- When the plasma (5 is large, MHD fluctuations are generally 
close to incompressible because pressure waves can travel rapidly enough to smooth out density 
disturbances. However, for short wavelengths, the radiation field decouples from the fluid, and 
Cs ^ Cg <^ VA, which means the plasma (3 becomes effectively quite small. For these short 
wavelengths, then, we can expect the turbulence to exhibit much greater compressibility. In the 
compressible regime, the speeds of the magnetosonic and Alfven waves are comparable, so they 
may couple much more easily. A similar effect happens for Alfven waves near recombination, as 
discussed by pubramanian and Barrow (1997)| . 

The slope and inertial range of the turbulent spectrum will also be affected by the plasma 
P parameter, which is usually held fixed in compressive MHD simulations ( [Matthaeus et ah] 



1996 ). Analytic theory and simulations show that for compressible MHD, 6p/p ~ {6va/cs) 



2 
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(where 5va = |(5B|/^47r/9), so when Cg drops dramatically in the non-diffusive regime, compressive 
damping will become very effective. Simulations of turbulent cascades with small (3 but 
with incompresible stirring will show how much energy can be transferred to compressible 
modes. Current simulations of compressible turbulence in the ISM (Charles Gammie, private 
communication) show that shocks form when ^ %, so that if the incompressible cascade does 
not transfer energy to compressible modes before reaching the non-diffusive scale, the energy 
may be dissipated in shocks at that scale. The dissipation in these shocks may be partly due to 
ordinary plasma processes, and partly due to radiation scattering. Thus, we expect that kmax will 
never be much greater than ko- As the radiation pressure varies with disk radius, fc/j changes and 
thus kmax changes, so the value of a may become a function of radius. 

Although certain consequences of radiation damping are relatively clear (at least qualitatively), 
consideration of this process also raises a number of questions: 

1) What is the nature of the coupling between compressive and incompressive modes? Is it 
large enough to allow the radiation damping rate to compete with the nonlinear frequency? Are 
the analytic estimates we have made useful in the nonlinear regime? 

2) In the simulations done to date, in which radiation pressure and transport are equally 
ignored, the magnetic energy density is an interesting fraction of the pressure and the associated 
fluctuations lead to a stress which is also proportional to the pressure. The question naturally 
arises whether, in radiation pressure-dominated disks, the r — (p stress and the energy in the 
magnetic field scale with the total pressure, or just with the gas pressure. The photon bubble 



instability ( Arons 1992 ) will likely affect the disk structure and stress (Gammie 1998). With 
explicit consideration of the quality of dynamical coupling between radiation fluctuations and fluid 
fluctuations, as outlined here, simulations should now be able to answer these questions. 

3) Can thermal or viscous instabilities be suppressed by radiation damping? Or does the 
dependence of dissipation on the radiation pressure exacerbate these instabilities? In both cases, 
the most important modes have radial wavenumbers < so the calculation here does not 
directly bear on them. However, one might expect that some of the same effects will qualitatively 
carry over. 

4) The relativistic portions of accretion disks may trap a number of long wavelength (i.e. 
kh < 1) normal modes (Nowak & Wagoner 1991, 1992). Some of these grow in amplitude due 
to viscous dissipation (Nowak & Wagoner 1992). Modulo the caveat of point 3), will radiation 
damping enhance (or destroy) these modes? 

5) Many seek the origin of disk coronal heating in the dissipation of rising MHD waves (e.g. 
Rosner, Tucker, & Vaiana 1978; Heyvaerts &: Priest 1989; Tout & Pringle 1996). If radiation 
damping quenches short wavelength fluctuations, will this affect the rate at which magnetic flux 
rises to the disk surface? 
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5.2. Gas-pressure dominated disks 

When gas pressure dominates over radiation pressure, radiation damping does not compete 
with the nonUnear frequency. The question of what causes the heating of the disk therefore 
remains open. This conclusion is equally true of conventional gas pressure-dominated disks and 
unconventional ones like ADAFs. 

Finally, the contrast between the radiation pressure-dominated and gas pressure-dominated 
regimes may mean that interesting observable effects occur in disks whose accretion rate fluctuates 
around the critical value of equation (^). If the value of a and the radiative efficiency depend on 
whether radiation damping plays a role, there could be significant modulations in the luminosity 
and spectrum on a viscous timescale. 
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